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CLAIMS 

What is claimed is: 

1 1. A method comprising: 

2 generating a first covariance matrix from a desired mean vector and a desired 

3 covariance matrix of a Bernoulli distribution; 

4 constructing a normal vector using the desired mean vector and the first covariance 

5 matrix; and 

6 generating a sampling vector using the normal vector and a threshold vector, the 

7 sampling vector having the desired mean vector and the desired covariance matrix. 

1 2. The method of claim 1 wherein generating the first covariance matrix 

2 comprises: 

3 generating an integral expression F for a first non-diagonal element Sy of the first 

4 covariance matrix at a row index i and a column index j, the integral expression having an 

5 integral limit as function of threshold elements Ti and xj in the threshold vector at the vector 

6 indices i and j; and 

7 obtaining the first non-diagonal element Sy using the integral expression F, a mean 

8 |o, k of the desired mean vector, and a desired non-diagonal element Sy of the desired 

9 covariance matrix. 

1 3. The method of claim 2 further comprising: 

2 obtaining a diagonal element sjj of the first covariance matrix at a first row index j 

3 and a first column index j using the mean uj at the vector index j, the diagonal element 

4 being equal to a desired diagonal element Ejj of the desired covariance matrix. 

1 4. The method of claim 3 further comprising: 

2 generating a threshold element xj of the threshold vector at a vector index j , the 

3 threshold element being equal to uj + Oj 4l inverf(l-2uj) wherein Uj and aj are desired 

4 mean and variance, respectively, at the vector index j and inverf is an inverse error 

5 function. 
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1 5. The method of claim 2 wherein constructing the normal vector comprises: 

2 generating normal elements of the normal vector using the desired mean vector and 

3 the first covariance matrix. 

1 6. The method of claim 5 wherein generating the sampling vector comprises: 

2 comparing a normal element Yk of the normal vector at a vector index k with a 

3 corresponding threshold element Xk of the threshold vector at the vector index k; 

4 setting a sampling element of the sampling vector at the vector index k to a first 

5 value if the normal element Yk is greater than the corresponding threshold element x^ and 

6 setting the sampling element of the sampling vector at the vector index k to a 

7 second value if the normal element Yk is equal to or less than the corresponding threshold 

8 element x k . 

1 7. The method of claim 2 wherein generating the integral expression F 

2 comprises: 

3 forming a first variable p = Sij/(GiO"j); 

4 forming a second variable c =-^2(l-p 2 ) ; 

5 forming a third variable P = (a, + aj)/(c yfl ), P being one of the integral limits; 

6 forming a fourth variable Q = (aj - aj)/(c -Jl ); and 

7 forming the integral expression 

8 F(p) = 4^V p2(1_P> (ofyfi+PiQ + P + P)~ erf(JU^(Q -p + P)))dp 

9 wherein p is an integral variable, erf is an error function, a; and aj are respectively 

10 equal to (T,-fij)/o-; and (Xj-u.j)/oj, i\ and Tj being the threshold elements at the vector indices 

11 equal respectively to the row index i and the column index j, p:; and p.j being the means at 

12 the vector indices equal respectively to the row index i and the column index j, <5\ and Cj 

13 being the variances at the vector indices equal respectively to the row index i and the 

14 column index j . 
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1 8 . The method of claim 7 wherein obtaining the first non-diagonal element 

2 comprises: 

3 determining a right hand side (RHS) quantity gy = + j-ii|aj; 

4 equating the integral expression to the RHS quantity to form an integral equation F 

5 = gij; and 

6 solving the integral equation for the first variable p. 

1 9. The method of claim 8 wherein solving the integral equation comprises : 

2 solving the integral equation using a numerical method. 

1 1 0. The method of claim 6 wherein the first value is 1 and the second value is 0. 

1 1 1 . A computer program product comprising: 

2 a machine useable medium having program code embedded therein, the program 

3 code comprising: 

4 computer readable program code to generate a first covariance matrix from 

5 a desired mean vector and a desired covariance matrix of a Bernoulli distribution; 

6 computer readable program code to construct a normal vector using the 

7 desired mean vector and the first covariance matrix; and 

g computer readable program code to generate a sampling vector using the 

9 normal vector and a threshold vector, the sampling vector having the desired mean 

10 vector and the desired covariance matrix. 

1 12. The computer program product of claim 1 1 wherein the computer readable 

2 program code to generate the first covariance matrix comprises: 

3 computer readable program code to generate an integral expression F for a first 

4 non-diagonal element sy of the first covariance matrix at a row index i and a column index 

5 j, the integral expression having an integral limit as function of threshold elements x { and tj 

6 in the threshold vector at the vector indices i and j ; and 
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7 computer readable program code to obtain the first non-diagonal element Sy using 

8 the integral expression F, a mean p. k of the desired mean vector, and a desired non-diagonal 

9 element Ey of the desired covariance matrix. 

1 13. The computer program product of claim 1 2 further comprising: 

2 computer readable program code to obtain a diagonal element Sjj of the first 

3 covariance matrix at a first row index j and a first column index j using the mean uj at the 

4 vector index j, the diagonal element being equal to a desired diagonal element Ljj of the 

5 desired covariance matrix. 

Q 

g 1 14. The computer program product of claim 13 further comprising: 

Si 2 computer readable program code to generate a threshold element Xj of the threshold 

W 3 vector at a vector index j , the threshold element being equal to Uj + Cj V2 inverf( 1 -2uj) 

4* ■ 

y 4 wherein Uj and Oj are desired mean and variance, respectively, at the vector index j and 

1. 5 inverf is an inverse error function. 

■O 
PJ 

ft 1 15. The computer program product of claim 12 wherein the computer readable 

y 2 program code to construct the normal vector comprises: 

l y 3 computer readable program code to generate normal elements of the normal vector 

4 using the desired mean vector and the first covariance matrix. 

1 1 6. The computer program product of claim 1 5 wherein the computer readable 

2 program code to generate the sampling vector comprises: 

3 computer readable program code to compare a normal element Y k of the normal 

4 vector at a vector index k with a corresponding threshold element x k of the threshold vector 

5 at the vector index k; 

6 computer readable program code to set a sampling element of the sampling vector 

7 at the vector index k to a first value if the normal element Y k is greater than the 

8 corresponding threshold element x k ; and 

9 computer readable program code to set the sampling element of the sampling vector 

10 at the vector index k to a second value if the normal element Y k is equal to or less than the 

11 corresponding threshold element x k . 



20 



04524.P061 

1 17. The computer program product of claim 12 wherein the computer readable 

2 program code to generate the integral expression F comprises: 

3 computer readable program code to form a first variable p = Sij/(o-jOj); 

4 computer readable program code to form a second variable c = V2(l-p 2 ); 

5 computer readable program code to form a third variable P = (ai + aj)/(c 4l ), P 

6 being one of the integral limits; 

7 computer readable program code to form a fourth variable Q = (aj - ai)/(c V2 ); and 

8 computer readable program code to form the integral expression 

9 F(p) = \e~ p2{1 - p) (erfJTT^iQ + p + P)- erf{^Tp{Q -p + P)))dp 

10 wherein p is an integral variable, erf is an error function, a; and aj are respectively 



11 equal to ("Ci-pi)/o-j and (T r Uj)/o-j, x\ and Xj being the threshold elements at the vector indices 

12 equal respectively to the row index i and the column index j, p* and uj being the means at 

13 the vector indices equal respectively to the row index i and the column index j, o"i and aj 

14 being the variances at the vector indices equal respectively to the row index i and the 

15 column index j. 



1 18. The computer program product of claim 1 7 wherein the computer readable 

2 program code to obtain the first non-diagonal element comprises: 

3 computer readable program code to determine a right hand side (RHS) quantity gjj 

4 = Zij + to; 

5 computer readable program code to equate the integral expression to the RHS 

6 quantity to form an integral equation F = gy; and 

7 computer readable program code to solve the integral equation for the first variable 

8 p. 

1 19. The computer program product of claim 1 8 wherein the computer readable 

2 program code to solve the integral equation comprises: 

3 computer readable program code to solve the integral equation using a numerical 

4 method. 

21 



04524.P061 



1 20. The method of claim 16 wherein the first value is 1 and the second value is 

2 0. 

1 21. A simulator comprises: 

2 a network modeler to model a network of free-space optical links; 

3 a reliability modeler coupled to the network modeler to evaluate a reliability model 

4 for the network; and 

5 a random sampler coupled to the network modeler and the reliability modeler to 

6 generate random samples for a Bernoulli distribution, the random sampler comprising: 

7 a covariance generator to generate a first covariance matrix from a desired 

8 mean vector and a desired covariance matrix of the Bernoulli distribution, 

9 a normal vector generator coupled to the covariance generator to construct a 

10 normal vector using the desired mean vector and the first covariance matrix, and 

11 a thresholder coupled to the covariance generator and the normal vector 

12 generator to generate a sampling vector using the normal vector and a threshold 

13 vector, the sampling vector having the desired mean vector and the desired 

14 covariance matrix. 

1 22. The simulator of claim 2 1 wherein the covariance generator comprises: 

2 an integral expression generator to generate an integral expression F for a first non- 

3 diagonal element Sjj of the first covariance matrix at a row index i and a column index j, the 

4 integral expression having an integral limit as function of threshold elements Tj and Tj in 

5 the threshold vector at the vector indices i and j; and 

6 a non-diagonal element generator coupled to the integral expression generator to 

7 obtain the first non-diagonal element sy using the integral expression F, a mean ja k of the 

8 desired mean vector, and a desired non-diagonal element Ey of the desired covariance 

9 matrix. 

1 23 . The simulator of claim 22 wherein the random sampler further comprises: 

2 a diagonal element generator to obtain a diagonal element Sjj of the first covariance 

3 matrix at a first row index j and a first column index j using the mean uj at the vector index 
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4 j, the diagonal element being equal to a desired diagonal element Sjj of the desired 

5 covariance matrix. 

1 24. The simulator of claim 23 wherein the random sampler further comprises: 

2 a threshold vector calculator coupled to the first normal vector generator to 

3 generate a threshold element x, of the threshold vector at a vector index j , the threshold 

4 element being equal to Uj + ctj a/2 inverf(l -2uj) wherein uj and cj are the desired mean and 

5 variance, respectively, at the vector index j and inverf is an inverse error function. 

1 25. The simulator of claim 22 wherein the normal vector generator generates 

2 normal elements of the normal vector using the desired mean vector and the first 

3 covariance matrix. 

1 26. The simulator of claim 25 wherein thresholder comprises: 

2 a comparator to compare a normal element Y k of the normal vector at a vector 

3 index k with a corresponding threshold element x k of the threshold vector at the vector 

4 index k; and 

5 a selector coupled to the comparator to set a sampling element of the sampling 

6 vector at the vector index k to a first value if the normal element Y k is greater than the 

7 corresponding threshold element tk and to set the sampling element of the sampling vector 

8 at the vector index k to a second value if the normal element Yk is equal to or less than the 

9 corresponding threshold element x k . 

1 27. The simulator of claim 22 wherein the integral expression generator 

2 generates the integral expression 

3 F(p)= ^-P °f e -P 2 d-p) (erffi+~^(Q + p + P)- erf(Jl7p~(Q -p + P)))dp 

4 wherein: 

5 p = Sij/((7i<Tj), c = a/2(1-p 2 ) , P = (aj + aj)/(c a/2 ), P being one of the integral limits, 

6 Q = (aj - ai)/(c V2 ), p is an integral variable, erf is an error function, a* and aj are 
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7 respectively equal to fa-pd/ai and (xj-jj.j)/crj, x { and xj being the threshold elements at the 

8 vector indices equal respectively to the row index i and the column index j, m and uj being 

9 the means at the vector indices equal respectively to the row index i and the column index 

10 j, <Ti and ctj being the variances at the vector indices equal respectively to the row index i 

11 and the column index j . 



1 28. The simulator of claim 27 wherein the non-diagonal element generator 

2 comprises: 

3 a right hand side (RHS) generator to determines a right hand side (RHS) quantity gy 

4 = Eij + mm; 

5 an equation solver coupled to the integral expression generator and the RHS 

6 generator to equate the integral expression to the RHS quantity to form an integral 

7 equation F = gy, and to solve the integral equation for the first variable p. 

1 29. The simulator of claim 2 8 wherein the equation solver solves the integral 

2 equation using a numerical method. 

1 30. The simulator of claim 26 wherein the first value is 1 and the second value 

2 isO. 

1 31. A system comprises: 

2 a processor; and 

3 a memory coupled to the processor, the memory having program code, the program 

4 code when executed by the processor causing the processor to: 

5 generate a first covariance matrix from a desired mean vector and a desired 

6 covariance matrix of a Bernoulli distribution, 

7 construct a normal vector using the desired mean vector and the first 

8 covariance matrix, and 

9 generate a sampling vector using the normal vector and a threshold vector, 

10 the sampling vector having the desired mean vector and the desired covariance 

11 matrix. 
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1 32. The system of claim 3 1 wherein the program code causing the processor to 

2 generate the first covariance matrix causes the processor to: 

3 generate an integral expression F for a first non-diagonal element Sy of the first 

4 covariance matrix at a row index i and a column index j, the integral expression having an 

5 integral limit as function of threshold elements x; and xj in the threshold vector at the vector 

6 indices i and j; and 

7 obtain the first non-diagonal element Sy using the integral expression F, a mean u. k 

8 of the desired mean vector, and a desired non-diagonal element Zy of the desired 

9 covariance matrix. 

1 33 . The system of claim 32 wherein the program code, when executed, further 

2 causing the processor to: 

3 obtain a diagonal element sy of the first covariance matrix at a first row index j and 

4 a first column index j using the mean Uj at the vector index j, the diagonal element being 

5 equal to a desired diagonal element Zjj of the desired covariance matrix. 

1 34. The system of claim 33 wherein the program code, when executed, further 

2 causing the processor to: 

3 generate a threshold element xj of the threshold vector at a vector index j , the 

4 threshold element being equal to p, + aj -J2 mverf(l-2uj) wherein uj and ctj are desired 

5 mean and variance, respectively, at the vector index j and inverf is an inverse error 

6 function. 

1 35. The system of claim 32 wherein the program code causing the processor to 

2 construct the normal vector causes the processor to: 

3 generate normal elements of the normal vector using the desired mean vector and 

4 the first covariance matrix. 

1 36. The system of claim 35 wherein the program code causing the processor to 

2 generate the sampling vector causes the processor to: 
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3 compare a normal element Y k of the normal vector at a vector index k with a 

4 corresponding threshold element x k of the threshold vector at the vector index k; 

5 set a sampling element of the sampling vector at the vector index k to a first value 

6 if the normal element Y k is greater than the corresponding threshold element T k ; and 

7 set the sampling element of the sampling vector at the vector index k to a second 

8 value if the normal element Y k is equal to or less than the corresponding threshold element 

9 Xk. 

1 37. The system of claim 32 wherein the program code causing the processor to 

2 generate the integral expression F causes the processor to: 

3 form a first variable p = Sij/(CTiOj); 

4 form a second variable c = ^2(1 - p 2 )); 

5 form a third variable P = (a; + aj)/(c V2 ), P being one of the integral limits; 

6 form a fourth variable Q = (a, - ai)/(c 4l ); and 

7 form the integral expression 

8 F(p) = ^—^[e-^-* (erfJlT^iQ + p + P)- erf(Ju r p~(Q -p + Ffftdp 

9 wherein p is an integral variable, erf is an error function, a s and a, are respectively 

10 equal to (Ti-Ui)/cTi and (Tj-Uj)/aj, n and tj being the threshold elements at the vector indices 

11 equal respectively to the row index i and the column index j, \x\ and uj being the means at 

12 the vector indices equal respectively to the row index i and the column index j, o"i and o-j 

13 being the variances at the vector indices equal respectively to the row index i and the 

14 column index j. 



1 38. The system of claim 32 wherein the program code causing the processor to 

2 obtain the first non-diagonal element causes the processor to: 

3 determine a right hand side (RHS) quantity gy = Sy + u^; 

4 equate the integral expression to the RHS quantity to form an integral equation F = 

5 gij; and 

6 solve the integral equation for the first variable p. 
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39. The system of claim 38 wherein the program code causing the processor to 
solve the integral equation causes the processor to: 

solve the integral equation using a numerical method. 



0. 



40. 



The system of claim 36 wherein the first value is 1 and the second value is 



